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ABSTRACT 

We study the emission of the short/hard GRB 090510 at energies > 0.1 GeV as observed by the 
Large Area Telescope (LAT) onboard the Fermi satellite. The GeV flux rises in time as t 2 and decays 
as t~ 1 ' 5 up to 200 s. The peak of the high energy flux is delayed by 0.2 s with respect to the main ~MeV 
pulse detected by the Fermi Gamma Burst Monitor (GBM). Its energy spectrum is consistent with 
F(y) oc v^ 1 . The time behavior and the spectrum of the high energy LAT flux are strong evidences of 
an afterglow origin. We then interpret it as synchrotron radiation produced by the forward shock of 
a fireball having a bulk Lorentz factor T ~ 2000. The afterglow peak time is independent of energy in 
the 0.1-30 GeV range and coincides with the arrival time of the highest energy photon (~ 30 GeV). 
Since the flux detected by the GBM and the LAT have different origins, the delay between these 
two components is not entirely due to possible violation of the Lorentz invariance. It is the LAT 
component by itself that allows to set a stringent lower limit on the quantum-gravity mass of 4.7 
times the Planck mass. 

Subject headings: Gamma rays: bursts — Radiation mechanisms: non-thermal — X-rays: general 



1. INTRODUCTION 

GRB 090510 is a short/hard burst at redshift 
2=0.903±0.003 (Rau et al. 2009) detected by Fermi 
(Guiriec et al. 2009), AGILE (Longo et al. 2009), Swift 
(Hoversten et al. 2009), Konus-Wind (Golenetskii et al. 
2009) and Suzaku (Ohmori et al. 2009). 

The Fermi-Gamma Burst Monitor (GBM) triggered 
on a precursor while the main emission episode in the 8 
keV-40 MeV energy range starts ~0.5 s after trigger and 
lasts up to ~1 s. The emission observed by the Fermi- 
Large Area Telescope (LAT) starts 0.65 s after the trigger 
and lasts ~ 200 s. The joint GBM-LAT spectral analysis 
showed the presence of two components. Fermi-LAT 
detected a 31±3 GeV photon delayed by 0.829 s with 
respect to the trigger (Abdo et al. 2009 - A09 hereafter). 

The precursor was not seen by AGILE, that triggered 
on the main emission episode. The flux detected by the 
Mini Calorimeter (MCAL, 0.3-10 MeV) lasts 0.2 s. As it 
ends, the Gamma Ray Imaging Detector (GRID, 0.03 - 
30 GeV) starts to detect a high energy component lasting 
for 10 s and decaying as (Giuliani et al. 2009). 

The distinguishing property of GRB 090510 arising 
from the Fermi and AGILE data is that the ~MeV emis- 
sion component, commonly detected in GRBs, is followed 
by a much longer lasting high energy emission detected 
above 100 MeV. Both the AGILE and Fermi spectra sug- 
gest that this component is not the extrapolation of the 
soft ^MeV spectrum to the GeV range. A09 interpret 
the ^MeV flux as synchrotron radiation and the LAT 
flux as its synchrotron self-Compton emission. The de- 
tection by Fermi of a 30 GeV photon sets a lower limit on 
the bulk Lorentz factor of the fireball T > 1000, based on 
the compactness argument (A09). The 30 GeV photon 
arrives 0.829 s after the trigger (set by the precursor) and 
0.3 s after the beginning of the GBM main pulse. These 
delays allowed A09 to put limits on the violation of the 
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Lorentz invariance. 

In this paper we propose a different interpretation of 
the emission detected by the Fermi-LAT and AGILE. 
If L > 1000 the fireball should start to decelerate and 
produce a luminous afterglow rather early (e.g. Piran 
2005), even at the sub-second timescale. By analyz- 
ing the Fermi-LAT light-curve and spectra we present 
strong evidences that the flux detected by the LAT is 
afterglow emission of the forward external shock. In this 
framework we derive the initial T of the fireball and set 
a lower limit on the quantum-gravity mass. 

Recent works on the high energy emission of LAT- 
detected GRBs include Kumar & Barniol Duran (2009), 
discussing GRB 080916C (Abdo et al. 2009a). Also for 
this burst they proposed that the LAT-detected flux can 
be synchrotron produced in the external shock (see also 
Gao et al. 2009 for GRB 090510; see also Fan et al. 2008; 
Zou, Fan & Piran 2009; Zhang & Peer 2009 for an inverse 
Compton origin). Hadronic models have been proposed 
to explain the emission properties of GRB 080916C (Raz- 
zaque, Dermer & Finke 2009). 

2. Fermi-L AT DATA ANALYSIS 

We have analyzed the Fermz-LAT data of GRB 090510 
with the Fermi ScienceTools (v9rl5p2) released on 
Aug. 8th 2008. Photons were selected (with the 
gtselect tool) around RA=333.552° and Dec=-26.598°. 
Different energy bins were considered for the analysis of 
the LAT light curve but only photons with energy >100 
MeV were extracted. Light curves and spectra were cre- 
ated with the gtbin tool. The spectral response files were 
created with the gtrspgen. The spectra were analyzed 
with Xspec(v. 12). 

3. RESULTS 

Fig. [T] shows the light curve considering all the LAT 
photons with energies >100 MeV. The times are scaled 
to T* = 0.6 s which corresponds to the time of the first 
main pulse observed by the GBM. We fit the light curve 
with the sum of two components, i.e. a smoothly broken 
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Fig. 1. — Ferrni-LAT light curve of the emission of GRB 090510 above 100 MeV. The times were scaled to the time T* = 0.6 s after 
the GBM trigger. This is the time corresponding to the main pulse of emission detected by the GBM in the 8 keV-10 MeV energy range 
(see A09). The solid line is the best fit to the data points obtained with a smoothly broken power-law plus a constant (dashed lines). The 
inset shows the AGILE light curve (photons energies between 30 and 300 MeV). The curve is the best fit of the Fermi data scaled to the 
AGILE data points. 

power-law and a constant to account for the flattening 
of the flux visible after 200 s: 

A (t/t h ) a 



R(t) = 



B 



(1) 



1 + (t/t h )<*+0 

When a > and (3 > Eq. Q] has a peak at £ pca k = 
i b (a/l3) 1 / {a+f3 \ The standard afterglow theory (Sari & 
Piran 1999) requires a = 2. Fixing a = 2, the best 
fit parameters (x 2 =14.6/14) are A = 385^0 counts/s, 
t pk = 0.217 ± 0.015 s, $ = 1.461^1 and B — A x 10~ 2 
counts/s. The best fit is shown by the solid line in Fig. 
[TJ The AGILE light curve (adapted from Giuliani et al. 
2009) of the photons detected by the GRID between 30 
MeV and 300 MeV is also shown in Fig. [1] (inset). The 
Fermi and AGILE light curves are consistent with the 
same decay law, i.e. i -1-5 , although AGILE missed the 
rising phase of the GeV emission. 

The emission above 100 MeV peaks at T - T*=0.22 
s (i.e. 0.82 s after the GBM trigger). The time of the 
peak coincides with the arrival time of the highest energy 
photon of 30 GeV. Fig. [JJ shows that the LAT flux lasts 
for about 200 s (and it sets to the background level after- 
words). Instead, the emission detected by the GBM in 
the 8 keV-10 MeV energy range ceases after ~1 s (A09). 

Fig. [2] shows the LAT light curve in the first 10 s 
separated in two energy bands, i.e. 0.1-1 GeV and >1 
GeV (top and middle panels, respectively). The curves 
correspond to the same best fit obtained from the >0.1 
GeV light curve (Fig. []}, only re-normalized to the data 
points. We further separated the light curve into four 
broad energy channels: 0.1-0.2 GeV, 0.2-0.4 GeV, 0.4- 
0.8 GeV and >0.8 GeV (Fig. [3]). Fig. [2]and Fig. [3]show 
that the time of the peak is the same in different energy 
ranges. 
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Fig. 2.— Eermi-LAT light curve of GRB 090510 between 0.1 
and 1 GeV and above 1 GeV (top and middle panels, respectively) 
in the first 10 seconds. The times are scaled to T*=0.6 s (sec 
text). The solid line is the fit of the light curve >0.1 GeV (Fig. 
[TJ- The bottom panel shows the photon spectral index (ler errors 
are shown) of the LAT spectra for the time-integrated spectrum 
(hatched region) and for three time resolved spectra (squares). 

We also analyzed the spectra of the early GeV emis- 



The GeV afterglow of GRB090510 



3 




0.1-0.2 GeV 



0.2-0.4 GeV 



0.4-0.8 GeV 



> 0.8 GeV 



0.1 



10.0 



1.0 
T-T* [s] 

Fig. 3. — Fermi-LAT light curve of GRB 090510 in four energy 
channels (from top to bottom): 0.1-0.2 GeV, 0.2-0.4 GeV, 0.4-0.8 
GeV, >0.8 GeV. The curves are the best fit obtained from the LAT 
light curve (>0.1 GeV - Fig. [TJ rescaled to the single channel light 
curves. 

sion component. We considered the spectrum integrated 
in time between T — T*=Q.l and 7 s and we also ex- 
tracted three time resolved spectra distributed in this 
time interval. The photon spectral index of the fit with 
a single power law of the average spectrum (hatched re- 
gion) and of the time resolved spectra (filled squares) are 
shown in the bottom panel of Fig. [5] The spectrum be- 
fore the peak is hard with a photon index 1.87±0.2 and 
then it softens to 2.2±0.2. Both are consistent with the 
time-integrated spectrum. 

4. ESTIMATE OF THE INITIAL BULK LORENTZ FACTOR 

The derived peak time of the LAT received flux trans- 
lates into an estimate of the bulk Lorentz factor Tq at 
the start of the afterglow. The peak time of the after- 
glow bolometric luminosity occurs at a time of the order 
of the deceleration time. If the circumburst number den- 
sity n is homogeneous we have (e.g. Sari & Piran 1999): 
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where -Ek,iso is the isotropic kinetic energy of the fire- 
ball, estimated through the emitted energy of the prompt 
emission assuming an efficiency rj (-Ek,iso = E^^o/rf). We 
use Eq. [5] to estimate T . Setting E 1: - mo = 3.5 x 10 52 
erg (A09, excluding the LAT component), rj = 0.2 and 
£ poak = 0.2 s we derive F = 1.96 x 10 3 rT 1 ^. This value 
is not much larger than the lower limits derived by A09 
through the compactness argument and assuming that 
the LAT component belongs to the prompt emission. It 
is also rather insensitive to the (unknown) particle den- 
sity n. The distance from the central engine correspond- 
ing to the peak time is i? pC ak ~ 2ci P eakFo/(l + z) = 
2.4 x 10 16 ™- 1 / 4 cm. 

4.1. A synchrotron origin of the LAT emission 

Following standard arguments, the minimum electron 
energy of the injected electrons in the forward shock is 
7, n ~ e e Tm v /m c , while the magnetic field value is B ~ 



r(87reB^w p c 2 ) 1 / 2 . Electrons with 7 m emit an observed 
synchrotron frequency u m ~ 2r(4/3)ei?/(27rm e c)7 2 a /(l-|- 
z) ~ 10.6r 4 (ne B ) 1/2 e 2 /(l + z) MeV. This frequency is 
below the LAT energy range, but the injection of a power 
law distribution of electrons extending to 7 ma x ~ (10 2 - 
10 3 )7m ensures that the LAT flux can indeed have a syn- 
chrotron origin. The synchrotron self-Compton (SSC) 
spectrum extends to much higher frequencies (e.g. Fan 
et al. 2008, Corsi et al. 2009), but becomes important 
only above v mt c ~ 7^m, i.e. above the TeV energy 
range. Note the strong dependence of v m and v m> c on 
the bulk Lorentz factor: a synchrotron origin of a ~GeV 
afterglow is reasonable only for rather large To, while the 
SSC flux becomes more important for smaller Tq. 

This has a simple and important consequence. Bursts 
with To ~ 100 or smaller can produce high energy af- 
terglow radiation through the SSC mechanism but the 
onset time of their afterglows will be large, in turn im- 
plying, for the same emitted energy, a lower luminosity 
(see also Kumar & Barniol Duran 2009). They are then 
more difficult to detect. The best candidates for a LAT 
detection are therefore bursts with large To, because this 
ensures an early onset of the afterglow, implying large 
luminosities. 

5. TEST OF THE LORENTZ-IN VARIANCE VIOLATION 

The arrival time of the 30 GeV photon coincides with 
the peak of the afterglow emission. This is reasonable, 
because this is the time when we have the maximum 
probability to detect it (maximum flux and hard spec- 
trum). If we assume that the 30 GeV photon was indeed 
produced at the peak time, then the maximum possible 
time delay it can have is of the order of the width of the 
time bin of the peak (0.15 s). More conservatively, we 
can assume that the 30 GeV photon was produced right 
at the beginning of the afterglow and it arrives delayed 
by 0.22 s due to violation of the Lorentz invariance. 

The time delay At between the arrival time of a low 
and a high energy photon in the case of a linear depen- 
dence of the photon's propagation speed on its energy 
is 



At = 



AE 



l + z 



M QG c 2 j Ho^n m (l + z) 3 + n A 



dz (3) 



where AE is the difference between the low and high pho- 
ton energy and Mqq is the quantum-gravity mass (A09; 
Amelino-Camelia et al. 1998; Jacob & Piran 2008). 

For a delay of 0.15 s we derive Mqq > 6.7Mpi anc k 
while the more conservative limit (delay of 0.217 s) is 
M QG > 4.7M Planck (we used h = 0.71, n A = 0.73, M = 
0.27). These limits are consistent with those of A09, but 
excludes their lowest estimates. 

6. CONCLUSIONS 

The detection of an early high energy emission in the 
GeV range inevitably flags a large bulk Lorentz factor 
To', to avoid suppression of the GeV emission due to the 
7-7 — > e ± process if the high energy photons belong to 
the prompt phase, or to have an early peak flux time if 
the emission belongs to the afterglow phase. 

We have shown that the latter case is indeed favored 
(see also Gao et al. 2009), because we see the peak 
time of the emission in the LAT light-curve (as also 
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seen in other GRBs in the infrared-optical band, see 
e.g. Molinari et al. 2007 for GRB 060418 and GRB 
060607A). Furthermore, also the energy spectral index 
F(v) oc v~ x is very similar to the ones we see in the af- 
terglow phase. A large To, implying an early onset of 
the afterglow, means a large luminosity at the peak time 
(for equal emitted energy), and large typical frequencies. 
This makes synchrotron the most likely process for the 
LAT emission we see. 

GRBs with smaller Tq will have their prompt emis- 
sion less blue-shifted, and it would be more difficult for 
them to reach the LAT energy range during their prompt 
phase. Their afterglows can, through the SSC process, 
but their afterglow peak time is longer, and so their fluxes 

are fainter (as t~^ a]l oc T^ 3 if they emit, at the peak, the 
same amount of energy, see Eq. 2). A large To, instead, 
means a large blue-shift for the photons of the prompt, 
an early onset of the afterglow, implying more flux at the 
peak, and finally larger intrinsic afterglow frequencies, al- 
lowing even the synchrotron photons of the afterglow to 
reach the LAT energy range. Therefore GRBs with large 
Tq should be much more luminous in the LAT energy 
range than the other GRBs (see also Kumar & Barniol 
Duran 2009). 

The limits derived here on the quantum gravity mass 



scale are not very different from the ones derived by A09, 
but we could associate the GBM and the LAT fluxes 
to two different components. We can then argue that 
the high energy photons are generated at (slightly) later 
times than the photons detected by the GBM, and the 
delay of their arrival times is not entirely due to quantum 
gravity effects. Instead, since photons above 100 MeV 
belong to the same component they are the best tool to 
investigate quantum gravity effects. 

This suggests a recipe for a robust test on the Lorentz 
invariance violation, possible with very bright and short 
bursts detected at high energies. A short duration of 
the prompt ensures that the fireball has a relatively nar- 
row width, and in turn this should correspond to a well- 
defined afterglow peak. A bright flux ensures good pho- 
ton statistics, enabling to measure more accurately pos- 
sible delays as a function of photon energies. 
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